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Abstract . 


The  paper,' which  is  the  first  in  a  series  of  two,  presents  an  approach 
which  allows  us  to  derive  a  family  of  homogenization  approaches  and  assess 
the  accuracy  of  any  homogenization  in  the  relation  of  given  input  data. 


1.  Introduction. 


-  The  study  of  periodic  media  is  one  application  of  partial  differential 
equations  that  have  highly  oscillatory,  periodic  coefficients.  Essentialy, 
the  problem  is  to  solve  th^  elliptic  differential  equation  , 

n 


(1) 


p.q=i 


lx)  =  f(x) 


on  n  c  with  prescribed  boundary  conditions,  in  which  a^^  and  a^  are 
real -valued  2it-periodic  functions  and  h  is  a  positive  number  that  is  small 
in  comparison  with  the  diameter  of  the  domain  n. 

The  problem  is  to  get  the  solution  of  (1)  for  relatively  (to  what?)  small 
h.  There  is  large  available  mathematical  literature  which  addresses  the  beha¬ 
vior  of  the  solution  of  (1)  as  h— ^0.  We  mention  here  for  example  [3],  [8], 

[ 9 ]  and  survey  [ 20 ] . 

One  of  the  main  applications  of  differential  equations  of  the  type  (1)  is 
in  the  field  of  composite  materials.  Here  the  aim  is  to  replace  the  composite 
by  homogeneous  materials  with  the  bulk  material  properties.  For  various 
ELspects  we  refer  to  [1],  [2],  [11],  [12],  [18].  A  brief  history  is  given  in 
[2].  The  accuracy  of  such  replacement  depends,  of  course,  on  the  goals  of  the 
analysis.  Hence  many  approaches  are  used  in  applications.  The  most  obvious 
approach,  namely  to  use  asymptotic  analysis  for  h— »0,  is  not  always  applic¬ 
able  because  h  is  given  and  ceinnot  be  changed  amd  because  of  particular  aims 
of  the  analysis.  In  a  similar  vein,  when  numerically  solving  the  problem  (1), 
one  faces  essential  difficulties  of  how  to  represent  the  microstructure  of  the 
compxsslte  materials.  This  difficulty  falls  into  the  claiss  of  solution  of 
elliptic  equations  with  rough  coefficients.  For  various  aspjects  of  this 
problem,  we  refer  to  [6],  [7]. 


As  was  said  above,  various  approaches  ceui  be  and  are  used  for  solving 
(1),  approaches  which  often  give  very  different  results,  see  e.g.  [10],  In 
addition  some  of  these  approaches,  although  in  principle  well  described,  are 
leading  to  leu'ge  technical  difficulties  because  a  lot  of  symbol  manipulation 
is  needed.  The  use  of  symbolic  manipulation  on  computers  does  not  simplify 
these  difficulties  too  much. 

This  paper  presents  and  thoroughly  aunalyzes  an  approach  which  is  directed 
to  overcome  the  various  major  difficulties  mentioned  above: 

a)  It  allows  the  design  of  aui  entire  class  of  "homogenization"  formula¬ 
tions  and  Judge  the  accuracy  and  reliability  of  any  homogenization  approach. 

It  also  allows  the  specification  of  the  class  of  problems  (e.g.  loads)  for 
which  a  homogenization  approach  is  applicable.  In  addition,  it  leads  to  a 
hierau'chal  construction  of  the  homogenization  formulations, 

b)  The  implementation  is  completely  numerical,  and  allows  adaptive 
modeling  (selection  of  the  equations). 

We  will  address  here  only  the  problem  with  H  =  although  very  impor¬ 
tant  features  of  the  solution  occur  near  the  boundary  when  12  is  a  bounded 
domain.  These  problems  have  a  special  character  and  will  not  be  addressed 
here.  Some  comments  will  be  made  in  section  5. 

We  will  assume  that 

i)  Q  =  r'' 

ii)  a^Ix)  i  >  0 

iii)  f  €  L2(r") 

and  that  the  problem  is  elliptic  and  self-adjoint. 

The  restriction  of  our  analysis  to  a  single  differential  equation  is  of  a 
technical  character  only,  as  the  ideas  are  also  applicable  to  a  system  of 
equations,  which  would  arise  in  elasticity  problems,  for  exajnple.  The  main 
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idea  of  the  approach,  under  the  assumptions  stated  above,  is  based  on  the 
result  that  the  solution  u^  of  (1)  can  be  written  in  the  form 


h,  ,  1 

u  (x)  =  - 


(2n) 


n/2 


f(t)0(^,h,t)e^'^*^dt 


in  which  f  is  the  Foiorier  transform  of  f  and  ^(y,h,  t)  is  a  function, 
which  is  2ir-periodic  in  y  and  analytic  in  h  and  t,  and  which  solves  the 
differential  equation 

(3)  ^  ^[^ap^(y)^(^(y)e^^^‘y)j  +h2aQ(y)^(y)  =  h^ 

p.q=i  ^  ^ 

on  {y  e  ;  ly^l  <  n} .  Other  representations  of  u^  that  are  related  to 
(2),  but  developed  in  a  different  context,  can  be  found  in  [8]  sind  [16]. 

By  taking  various  expansions  of  ♦  with  respect  to  h  and  t,  we  can 
express  approximately  the  solution  u^  in  terms  of  solutions  of  auxiliary 
partial  differential  equations  with  constant  coefficients,  or  even  pseudodif¬ 
ferential  equation,  or  alternatively,  we  cam  design  a  system  of  "ansatz"  func¬ 
tions  to  be  used  in  a  finite  element  model.  Considering  the  error  estimates 
associated  with  the  expansions  of  0  allows  us  to  design  an  adaptive  method 
of  selecting  a  "model"  that  would  yield  an  approximate  solution,  whose  accu¬ 
racy  meeds  a  prescribed  tolerance.  These  ideeis  are  more  fully  discussed  in 
[6],  where  we  introduced  a  method  to  systematically  derive  numerical,  computer 
oriented  methods  for  an  approximation  of  u^. 

In  this  paper,  we  concentrate  our  attention  on  the  representation  (2), 
whereas  we  will  m2dce  a  thorough  auialysis  of  ^(y,  h,  t)  in  [15].  Consequently, 
the  properties  of  0  that  are  used  in  this  paper  will  be  stated  without 
proof.  The  integral  in  (2)  is  defined  8is  a  Bochner  integral  of 

-valued  function  is  defined  in' the  next  section).  As  a 
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simple  application  of  (2),  we  will  given  an  alternate  proof,  in  section  5,  of 
the  classical  homogenization  result  (the  limit  of  as  h — >0). 

This  paper  and  [15]  are  based  on  the  first  author’s  Ph. D.  thesis  [14],  in 
which  additional  details  and  references  csm  be  found. 


2.  Notation  and  Statement  of  the  Problem. 

For  J  =  0, 1,  and  for  any  w  e  R,  define  the  weighted  Sobolev  space 

Hj|(IR*’)  to  be  the  completion  of  Cq(Ir”)  (the  complex-valued  C*-functions 

that  have  compact  support  on  R^),  with  respect  to  |1*||  .  ,  where 

0  >  ^ 

Z,  _,a  ,  V  ,  2  2w  I X 1  . 

ID  u(x) I  e  dx. 

R"  laUj 

(For  X  €  R*',  1x1  3  lx,  1  +  .  .  .  +  lx  1 .  )  We  will  use  H'^(r”)  and  11*11  ,  to 

In  J 

denote  the  standard  Sobolev  space  and  norm  on  R^  (i.e.,  when  y  =  0).  Next, 
we  Introduce  the  Sobolev  spaces  of  periodic  functions  for  which 

S  s  (y  s  (y^ . y^)  €  r”  :  lyj^l  <  n  for  k  =  l,...,n> 

is  the  fundamental  period.  For  J  =  0, 1,  we  denote  the  standau'd  Sobolev  norm 

on  S  by  11*11,  _•  we  define  (S)  to  be  the  completion,  with  respect 

J I b  per 

to  the  norm  11*11,  of  the  complex-valued  c"-f unctions  on  R^  that  are 
J  »  ^ 

2n-periodic  in  each  coordinate  variable. 

Let  a  (p,q  =  l,...,n)  and  a_  be  real-valued,  2Tt-periodic, 

Pm  v 

L  -functions  defined  on  r”.  Furthermore,  assume  a  =  a  and  aissume  that 

oo  qp  pq 

there  exist  positive  constants  and  such  that 


4 


(4) 


a^Cx)  >  tQ  and 


n 


n 


Z  a  (x)<<  ir  JKr  for  all  <  a  C, 

lp.q=l  ^  ^  ^  ^p=l  ^  ^ 


almost  everywhere  on  r”.  For  each  h  >  0,  define 

/-  n 


4'(h)[u,v]  = 


‘P.q-l  ”  <1  P 


klx. 


An  immediate  consequence  of  the  conditions  imposed  on  the  coefficients  a^ 
and  a^^  is  that  there  exists  a  constant  C,  independent  of  h  >  0,  such 


that 


(5) 


14'(h)[u,vl|  S  C||ull^||v||^  and 
H'(h)[v,v]|  a  min{yQ,yj>||v||j 


for  all  u  and  v  in  H  (R  ).  Then,  according  to  the  Laoc-Milgram  theorem, 
for  each  h  >  0  and  each  f  €  L^,  there  exists  a  unique  function  u^  € 
that  satisfies 


(6) 


♦(hltu  ,v]  = 


f (x}v(x)dx  for  all  v  e  H^CR*^), 


because 


(7) 


f (x)v(x)dx 


,n 


is  a  bounded  linear  functional  on  H^(R^) 


Next,  for  each  h  e  C  and  t  €  C*',  define  the  sequi linear  form 

4Kh,t)  :  (s)xH^  (s)-^C  by 

per  per 


♦(h,t)[^,v]  s 


r  n 


^  apq(y)^(0(y)e^^^’^)^(v(y)e  ^^'^’^}+h^aQ(y)^(y)v(y)  idy 


S'p,q=l 
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n+ 1  n+ 1 

Lemma  1.  A  neighborhood  G  c  C  ,  of  IR  cam  be  found  such  that  for  each 

(h,  t)  e  G,  there  exists  a  unique  fvinction  ^(‘.h,  t)  6  (S)  that  satisfies 

per 


«(h,t)[^(*.h,t).vl  =  h^ 


vTyTdy  for  all  v  €  (S). 


per 


Furthermore,  the  mapping  (h,t)  €  G  i — »^{‘,h, t)  €  (S)  is  holomorphic,  by 

which  we  mean  that  about  each  point  in  G,  the  function  (h,  t)i — >0(»,h,t) 

can  be  expeuided  in  a  power  series,  convergent  in  which  each 

coefficient  is  am  element  in  (S). 

pier 

For  the  most  pairt,  the  proofs  of  statements  concerning  ^(•,h,  t)  are 
omitted  in  this  paper  since  we  give  a  fairly  complete  analysis  of  ^(•,h,t) 
in  [15]. 

In  section  4,  we  show  that  u^  admits  the  representation 


(8) 


u  (x) 


(2n) 


n/2 


f(t)(>(^,h,  tle^^’^^dt, 


in  which  f(t)  =  - ^ ^  f  f(x)e  ^dx  and  in  which  the  integral  is  a 

(2n)^^  Jr'' 

Bochner  integral  of  -valued  functions.  Our  proof  of  (8)  ham  as  its 


first  step  the  claim  that  for  each  h  >  0  auid  t  e  r". 


(9)  XI— »^(^,h,t)e“‘’' 

n 

solves  (6)  when  f(x)  =  e^^  However,  (9)  is  not  am  element  of  H^(r''), 
aind  for  this  choice  of  f,  (7)  is  not  a  bounded  lineair  functional  on  H^Cr"). 
Consequently,  we  consider  ♦(h)  as  a  sesquilineam  form  on  H^^(r'')xH^(r'') 
for  (sufficiently  small)  pxjsitive  numbers  v. 

The  main  tool  for  ainalyzing  ’Kh)  is 

Theorem  2.  Let  and  be  two  complex  Hilbert  spaces  with  respective 
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norms  ll’llj^  and  associated  Inner  products  ^  Let  B[*,*] 

be  a  sesqulllneeir  form  defined  on  H^xH^  for  which  there  exist  positive  con¬ 
stants  M  and  y  such  that 

a)  |B[u,v]I  ^  Mllull ^||v||2  for  all  u  e  and  v  e 

b)  inf  sup  |B[u, v]|  >  y  >  0,  and 

ueHj 

l|u||j=l  Hull 2^1 

c)  sup  |B[u, v]l  >  0  for  each  v  6  H„,  v  *  0. 
ueH^ 

If  f  e  H2,  the  space  of  bounded  conjugate-lineau'  functionals  on  then 

there  exists  a  unique  Uq  €  such  that 

d)  B[Uq,v]  =  f(v)  for  all  v  e  H2.  and 

e)  IIUqIj  ailfl  .. 

A  proof  of  theorem  2  in  the  case  of  real  Hilbert  spaces  cam  be  found  in 
[4]  (as  Theorem  5.2.1).  The  method  of  proof  in  the  complex  case  is  essen¬ 
tially  unchanged  and  thus  will  be  omitted. 

We  now  prove 


Lemma  3. 

There  exist  positive  constants  ^q.C, 

and  y 

such  that 

V  e 

(O.Uq 

)  and  al  1 

h  >  0. 

i) 

I'Kh)  [u,  v] 

ii) 

inf 

sup  (4'(h)(u, vJI  >  y  >  0, 
llv||i^^=l 

and 

iii) 

sup  I'Kh)  [u,  v]  1  >  0  for  all  v  € 

ueH^^dR*') 

hJ(ir”) 

and  V  ^  0 

The 

constants 

and  y  aire  indepjendent  of 

h  >  0; 

however , 

on 

‘'o- 

Proof.  Statement  (i)  follows  because  4'(h)  has  L^-coeff icients.  To  prove 
statement  (ii),  define  (Tu)(x)  »  u(x)e  for  u  e  Then  for 
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u  >  0, 


|2  ,  [y,au_ 

'l.y  ^  Z-'Sx, 


(x)  ~2v  sgn(x  .)u(x)  1^  +  lu(x)  1^  e 


S  2(l+4nw^)||u||^  . 

i.  >  1/ 


4'(h)lu,Tu]  = 


r  r  ^ 

Z.X.3U 

^pq  h  3x 
11  a 

•’lR“4),q=l  ^ 


(x)-g^^  (x)  +  aQ(^)u(x)v(x)|dx 


IR  ^,q=l 


q  P 


-2y |x| 


=  -2v 


•  r  ” 
IR''^.q=l 


a  (~)|^(x)sgn(x  )u^x)  « 
pq  h  ox  p 

q  ^  J 


s  4'j(h)lu]  -  2i;4'2(*'nu]. 


A  simple  consequence  of  (4)  is  4'^(h)[u]  2  min{9rQ, ||u|| ^  We  also  have 

2 

S  c||u||j  ^  for  some  constant  c,  Independent  of  h.  Combining 

these  two  inequalities  with  (10)  zuid  (11)  yields 

min{rn»3r,>  -  2cv 

[♦(hllu.Tu]  I  >  (min(3f„,r,>  -  2cv)l|u||7  - Hull,  .JlTuH, 

U  1  1,-p  ! - ^  l,-v  l.y 

v2(l+4nu  ) 

Consequently,  there  exist  positive  constauits  amd  y,  which  are  indepen¬ 

dent  of  h  >  0,  for  which  't'(h)[u, Tu]  >  yl|u||^  _^l|Tu|l^  ^  for  all  v  e  (O.^q) 
This  proves  (ii). 

Statement  (iii)  is  proven  in  a  similsir  memner.  □ 


Throughout  the  remainder  of  this  paper,  we  implicitly  assume  v  e  (O.u^) 
Let  f  6  L^dR*^);  then 
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I  f(x)vTOdx|  S  ||f|LI|v||^  ^  for  all  v  e  hJ(IR"). 

•  r" 


Lemma  3,  in  conjxinction  with  Theorem  2,  now  yields 

Theorem  4.  For  each  h  >  0  and  f  €  L2(R^)  there  exists  a  unique  function 

u^  €  (r'^)  for  which 

~v 


(12) 


"KhlCu  ,v]  = 


f (x)v(x)dx  for  all 


V  €  H^(r’^). 
V 


Furthermore.  Hub'll ^  i 

There  is  no  ambiguity  in  denoting  the  unique  solutions  of  (B)  and  (12)  by 
u^  because  H^(R^)  c  H^(R^)  c  H^^(R^)  implies  they  are  the  same  function. 
However,  (12)  can  be  solved  when  f  belongs  to  a  broader  claiss  of  functions 
than  the  class  of  L2-functions;  namely,  when  f  belongs  to  the  dual  space  of 
One  such  function  is  defined  by  f(x)  =  e^^ 


3.  Prel  iminau’ies. 

It  is  easy  to  prove 

Theorem  5.  There  exists  a  constant  Cq  independent  of  t,  such  that 

llze^’'''y|ll,S  5  CQ(l+|lt||)||;tllj_g 

for  al  -  X  €  H^(S)  and  for  all  t  €  r”,  where  ||t||^  s  tf  +  .  .  .  +  t^. 

1  n 

For  each  h  >  0  and  each  w  €  z”,  define 


S(h,u)  a  (x  €  R^ 


(«j-l)nh  < 


(Wj+l)ith, 


J  =  1 . n}. 


Note  that  S( 1, (0, . . . ,0) )  =  S  is  the  fundaunental  domain  for  the  periodic 


spaces . 


Lemma  6.  For  each  h  >  0  aind  y  >  0,  there  exists  a  constant  C^^Ch.y)  that 
remains  bounded  eis  h— >0,  such  that  for  all  t  e  r”,  zuid  for  any  %  e 

1)  s  Cj(h,..)|zlo  s 

and  for  any  x  «  (S), 

per 


ii)  ^  (1+h  ^)Cj(h.y)llze^^’'‘^||j_5. 


Proof. 


M  it 

llz(jj)e 


•X,-.  ^  i„  1  d 

W€Z  *’s(h.2u)  >15=1 


,  8  ,  ,x.  it*x.,2^  I  ,x.,2  -2y|x|  . 

I^(z(j^)e  )|  *\xi^)\  e  dx. 


where  2w  a  (.2u. ,  .  .  .  ,2u>  ).  Making  the  substitution  ^  -  y  +  27rw  in  the 
in  h 

integral  over  S(h,2a))  auid  using  the  periodicity  of  x  yields  a  constant 
c,  independent  of  h,  such  that 


,,  ,x,  it*x,,2  V*  fu'lV'i  3  f  (  \  i2n:hwt,2^  ,  ,  .,2. 

”i.-w 2-  ‘  r  Z^’ay“  ^  ^  ®  ®  +lz(y)l  ) 


uez"^''s  p=l  P 


-2uh|y+27ruL  n . 
•e  ^  h  dy 


.  -..^.-2,  -4iwh|w|  .n,,  iht*y„2 

^  c(l+h  )  )  e  h  Wxe  g. 


It  is  not  difficult  to  prove  that  there  is  a  constemt  c  that  is  independent 
of  h,  such  that 


.n  -4iivh|u|  .  .  n, , 

h  /  e  i  ch  ( 1+  e 


11/ , 

1  ( 1+  e 


-rm/hlxl 


Upon  setting  Cj(h,u)  a  v  cc(h’^+(2ny)  ”),  (ii)  follows. 
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A  similEir  argument  in  which  the  contribution  from  the  first  order  deriva¬ 


tives  is  ignored,  produces  (i). 


Lemma  7.  For  all  h  >  0  2uid  t  e  we  have 


ll^(-.h.t)|L  o  ^ 


(27t) 


n/2 


ajid 


O.S  -  y. 


ll^(-.h.t)e^^^*yilj  g  i  CgCh) 


for  some  positive  number  CgCh). 


Proof.  It  follows  from  (4)  and  Lemma  1  that 


^^I^C^Cy.h.tle^^'^'yil^dy  +  yoh^ll^C-.h.t)!!^  , 


S  p=l 

5  *(h,t)[0(«.h,t),«(*.h,t)] 


S  (2ii)"^V||^(«.h.t)||p  g. 


The  lemma  now  follows  with  CgCh)  s 


(2n) 


n/2 


_2  ' 

min{yo.rih  } 


4.  The  representation  of  u^. 
We  begin  with 


Theorem  8. 

For  each  h  >  0 

auid  t 

ir.n 

€  R  , 

i) 

,fX  ,  .  %  it*x 
XI — h,  t)e 

is  in 

H^^(r"),  and 

ii) 

'*'(h)[0(^,h,t)e^’^*’‘, 

V, 

Jfl 

e^^  %(x)dx  for  all  v  €  H^(R^) 

jtt  V 

Furthermore, 
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yu 

Proof.  Statement  (i)  follows  from  Lemmas  6  eind  7,  which  imply 

h,  t)e^^  ^11  j  ^  (1+h  ^)Cj^(h, y)C2(h).  Assuming  that  (ii)  is  true,  (iii) 

is  a  consequence  of  Theorem  2,  Lemma  3,  suid  the  fact  if  e^^  ^v(x)dxl  6 

Jr” 

— ^||v|L  .  The  proof  of  (ii)  is  based  upon  determining  the  relationship 

n/c  1, V 

V 

between  ’tih)  auid  #(h,  t),  and  then  using  Lemma  1. 

00 

For  each  h  >  0,  there  exists  a  locally  finite,  C  -peu'tition  of  unity 
{<r  (‘.h)  :  u  e  Z^>  subordinated  to  {S(h, w)  :  w  6  Z^}  such  that 

Ct) 

r  <r  (•,h)v  converges  to  v  in  whenever  v  €  H^(IR^).  The  basic 

_n  u  V  V 

ueZ 

9<r 

requirement  of  the  partition  of  unity  is  that  |<r  (x,h)|  and  |^— (x,h)|  for 

(i>  ox 

P 

p  =  l,...,n  are  uniformly  bounded  for  x  6  and  u  6  Z^.  Then  v^(*,h)  s 
<r  (^hlv  ahs  compact  support  in  S(h, w),  auid  for  emy  x  ^  (S), 

Cl)  p©p 


'Kh)(y(g)e^’^*^,v]  =  ^  *(h)  v^(x,  h)  ] 


U€Z 


n  J 


/■  n 


Z,x,  d  ,  ,x,  it'X.^'^w,  .  , 

S(h,  w) '■p,q=l  ^ 


it  'x 


V  (x, h) 

U) 


dx 


1  n  In 

because  'Kh)  is  a  continuous  sesquilinesu*  form  on  (R  )xH^(IR“),  accord- 

-V  V 

ing  to  Lemma  3  (i).  In  an  effort  to  deform  the  region  of  integration  S(h,(j) 
into  S  =  S(  1,  (0, .  . . ,  0) )  in  each  integral,  meike  the  substitution  ^  =  y  +  irw 
for  X  6  S(h,u),  in  which  «  is  the  n-tuple  of  even  integers,  that  is, 
derived  from  u  according  to 
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(J  = 

p 


u> 


(j)  even, 
P 


.0)  -1,  w  odd. 

V  p  '  p 


Using  the  periodicity  of  a_,  {a  :  p, q  =  1 . n},  and  x>  each  integral 

u  pq 

in  (13)  becomes 


(14) 


S(l.w-u) 


n 


av. 


h'^  ^  ap^(y)^(z(y)e^^^*^^'*’""^)^(h(y+7ru),h) 

p.q=l  ^ 


+  aQ(y)3c(y)e^^*'*  ^^^”“v^(h(y+7ru) ,  h) 


i-Oj 

•h  dy. 


Next,  for  each  6)  €  Z  ,  define 


(15) 


v^(y, h,  t)  s  V  (h(y+Jrw),h)e  (y+ww) 

U  ti 


for  y  €  S(l,u-u)  and  extend  v^(»,h, t)  to  all  of  by  2ir-periodicity. 

u 

Since  the  support  of  yi — >v  (h(y+nw),h)  is  contained  in  S(l,w-u),  it 

u 

follows  that  v^(*,h,t)  e  (S).  Lfeing  (15)  to  substitute  for 
w  per 

v^(h(y+Ttu) , h)  in  (14)  yields 


S( 1 , w-w) 


n 


,-2  f  \  d  t  t  \  iht*y«  9,0,  ,  .  -iht*y, 

p.q*i  ^ 


Vh^dy. 


Now,  the  domain  of  integration  S(l,«-w)  can  be  replaced  with  S(1,(0 . 0)) 

=  S,  auid  consequently 


(16) 


4'(h)[z(r)e^^*^,v]  =  Y'  h”  ^♦(h,  t ) [}c,  v°( • , h,  t )  ] 

h  Z^n  " 


ucZ 


for  all  V  €  hJ(r”). 
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Noting  Lemma  1,  It  is  now  a  simple  matter  to  prove  (ii): 


’Kh)[0(J,h.t)e^^*^.v]  =  y  h'^"^«(h.t)[0(*,h,t),v°(-,h,t)] 
h  u) 


ueZ 


=  1/ 


(deZ 


V  (y.h,t)dy 

W 


it*xv(x)dx 


,^n 


for  all  V  e  H^(IR”),  since  v*^(*,h,t)  €  (S). 

V  u  per 


Lemma  9.  For  each  h  >  0,  t  e  IR*'i — t)e^^  ^  6  is  a  continuous 

mapping. 

Proof.  The  continuity  of  ti — )e^^^  e  follows  in  a  straightforward 

manner. 

Upon  setting  t  =  0  and  %  =  ^( • ,h, t )  -  ,  h,  t)  in  Lemma  6 

limll^(ph.t)  h.T)|| 

t-»T  ’ 

<  (l+h"^)C.(h.i;)liml|^(-,h.t) -^(•.h,T)|L  „ 

^  t-»T 


=  0 


follows  from  Lemma  1. 


For  each  f  €  L2(IR^),  the  Fourier  treuisform  of  f  is  defined  by 


f(t)  = 


(2rt) 


n/2 


.  -it*x, 
f(x)e  dx. 


The  following  notation  will  be  used  with  the  function  f  only  (as  it  appears 

in  (12)).  For  any  f  €  L2(IR'')  and  N  >  0.  define  f^  e  L2(IR^)  as  the 

inverse  Fourier  transform  of  fL.  ,  i.e. 

'it:  ||t||SNl 
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f„(x) 


(2ir) 


n/2 


A  i  fc  •  V 

f(t)e  ^  dt 


and 


ItllSN 


lltll  i  N 
lltll  >  N. 


Parseval's  inequality  implies 


limllf-f^llQ  =  limllf-f^llg  =  0. 
N-»fl9  N-»« 


Next,  for  each  N  >  0,  define  u^(»;N)  €  by 


u^(x;N)  =  -  ^ 


(2n) 


n/2 


f(t)^(g.h,t)e^^'’^dt 


ItlUN 


fj^(t)^(^  h,t)e^^*Vj, 


in  which  the  integral  is  to  be  interpreted  ais  a  Bochner  integral  of 
valued  function  (cf.  [17]).  In  order  to  show  that  the  integral  has  such  a 
meaining,  we  need  to  show  that  the  integrand  is  strongly  measurable  aaid  that 
ti — »||f(t)^(^,h,t)e^^  ^llj  is  Lebesgue  integrable  over  {t  ;  ||t|l  i  N} . 

This  integrabillty  condition  is  satisfies  ais  a  result  of  (iii)  of  Theorem  8 
and  because  f  is  an  L^^-f unction  on  {t  :  ||t|]  i  N}.  By  strong  measurability 
of  the  integrand,  we  mean  there  exists  a  sequence  of  simple  functions 
ti — >Wj^(*,t)  €  such  that 

(17)  liml|f(t)^(p,h.t)e^’^'’^- w  (x,t)l|  =0. 

k-»09  "  K  1,-v 


It  is  eaisy  to  construct  such  a  sequence  (and  we  do  so  in  the  proof  of  Theorem 
10  below)  using  Lemma  9  and  the  measurability  of  f.  Furthermore,  if  the 
sequence  of  simple  functions  satisfies 


(18) 


lim 
k-Ka  J 


||f(t)0(^,h.t)e^’^' 


-  W|^(x,t)| 


1.-1^ 


=  0 


ItlUN 
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then 


Proof.  The  essence  of  this  proof  is  that  a  sequence  of  simple  functions 

satisfying  (17)  and  (18)  can  be  chosen  so  that  each  value  of  each  simple  func- 

V  It  •  X 

tion  is  of  the  form  h,T)e  for  some  t  €  {t  :  ||t||  <  N}  aoid  for  some 

complex  number  c.  With  this  being  the  case,  Theorem  8  (ii)  caji  be  used  to 
evaluate  ♦(h) (u,v]  whenever  u  is  the  integral  of  one  of  these  simple 
functions. 

To  begin,  note  that  H^(Ir”)  c  Lj(IR^).  Consequently,  v  is  continuous 
(cf.  [15])  when  v  e  H^(IR^),  and  it  follows  from  Theorem  8  that 

V 

♦(h)[0(^,h,T)e^‘^*^,v]  =  (2n)'^^0(T) 
is  well  defined  for  each  t  e 

Set  n  H  (t  e  :  ||t||  £  N>.  Now  select  a  sequence  of  simple  functions 
Sj^  :  n— )C  and  a  sequence  of  collections  ^  :  J  =  0,...,Mj^>  of  measur¬ 

able  subsets  of  £1  for  k  =  1,2,...  and  for  which  the  following  is  true. 

For  each  k: 


a) 

JsrO  ’ 

“k,j 

b) 

for  J  =  1 , .  . 

.  ,  Mj^,  if  t  and  t 

are  in 

"k.j’ 

V  1  t  ^ 

ll^(^,h,t)e'’' 

’X  ,,x  .  .  .  it»x,, 

<  2"*', 

and  |v(t)  -  v(t) |  < 

2  and 
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c)  Sj^  is  constant  on  j  for  J  =  0 . with  =  0  for 

t  e  Qj^  Q,  Sj^— »fj^  pointwise  a. e.  on  n,  and  |Sj^(t)|  i  |fj^(t)| 
a.  e.  on  n. 

k  1 

Next,  for  each  k,  pick  a  t  e  ^  for  each  j  =  1 . Mj^,  and  define 

ro,  t  .  n, 

'  “k,J’ 


Wj^(x,t)  = 


.fX  .  k,  J.  ix'^'  '^•x 

^(jj.h.T  ’'')e 


and 


rk(t)  = 


"  Sc.O 

t  €  J. 


A  consequence  of  the  preceding  construction  eund  of  Theorem  8  (iii)  is 


lf^(t)^(^.h.t)e^*'*’'-Sj^(t)Wj^(x.t)||j 


7V 

+  2*‘'l|f„(t)|. 


It  follows  from  the  first  inequality  that 


lim!|fj^(t)0(g.h.t)e^’^*^-Sj^(t)Wj^(-,t)||j  =  0  a.e. 

k“^oo  * 


Furthermore,  the  Lebesgue  dominated  convergence  theorem  auid  the  second 
inequality  imply 


lim 

k-»» 


|fj^(t)0(^,h,t)e^^'^-Sj^(t)Wj^(.,t)||j  _^dt  =  0. 


R 


Consequently,  u^'C^N)  =  lim(2n)"‘^^r  s  (t)w  (•,t)dt  in  (IR"). 

k-»«  ■'lltllSN  ■*' 

Finally,  the  continuity  of  ui — >^(h)[u, v]  (Lemma  3);  the  definitions  of 


17 


and  and  Parseval’s  equality  implies 


'l'(h)tu^(‘:N),v]  =  lim(2it)"'^^ 

k->oo 


Sj^(t)Wk(*,t)dt,v 


=  lim 
k-»oo 


Sk(t)rk(t)dt 


|tll<N 


fj^(t)v(t)dt 


|t||<N 


f.,{x)v(x)dx 

N 


for  all  V  e  because  =  0  for  ||t||  >  N. 

V  N 


It  is  now  a  simple  matter  to  prove  the  main  result  of  this  paper. 

Theorem  11.  Suppose  h  >  0  and  f  €  LgdR^).  Let  u^  €  be  the 

solution  of  (12).  Then 


(19)  u^(x)  =  lim(2Ti)“^^ 

N-»oo 


f(t)^(^,h,t)e^^‘^dt  in  (R^), 
h  -r 


|tll<N 


where,  for  each  N,  the  integral  is  defined  as  a  Bochner  integral  of 
H^^(R^) -valued  functions. 

Proof.  A  consequence  of  Theorem  10  is  that 


(20) 


«'(h)[u^- u^(*;N),v]  = 


(f (x)-f„(x) )v(x)dx 
N 


for  all  V  e  H^(R),  from  which  the  inequality 


(21) 


is  easily  derived  (cv.  Theorem  2  and  Lemma  3).  Then  (19)  follows. 
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Actually,  (19)  converges  in  H^(IR*')  even  though  each  integral  is  defined 
only  as  a  function  in  The  reaisoning  that  allows  us  to  identify  the 

unique  solutions  of  (6)  and  (12)  also  yields  u^(*;N)  e  H^dR*^)  and  the  fact 
that  (20)  is  valid  for  all  v  e  H^(1R”).  Now,  the  Lax-Milgram  theorem  suid  (5) 
imply 


I  h  h, 
|u  -  u  ( 


^“>"1  ^  STHTV^"''- Vo- 


5.  Homogenization. 

In  this  section  we  derive  first  the  classical  result  of  homogenization 
which  states  that  u^  converges  as  h  tends  to  0,  to  a  function  that  is 
the  solution  of  a  constant-coefficient  partial  differential  equation.  This  is 
an  example  of  analyzing  u*'  through  (19)  emd  an  analysis  of  t). 

According  to  Lemma  1,  it  is  possible  to  expand  0(*,h,t)  in  powers  of 
h,  for  each  t  €  IR*^.  Consequently,  we  can  write 

(22)  ^(‘.h.t)  *  ^Q(*.t) +^j(-,t)h+ . .  .  . 


The  functions  {^j(*,t)  :  J  =  0,1,...}  can  be  determined  by  expanding  (3)  in 
powers  of  h  and  substituting  (22).  Here,  we  ane  interested  in  only  the 
constant  term,  and  solving  for  it  yields 


^o(-.t)  =  ggCt)  - 


y  A  t  t  +A^ 

p.^^t  w  'I  p  0 


where  A„  and  { A 

0  pq 


p.q 


l,...,n}  ane  derived  from  the  periodic  coeffi- 
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cients  and  {a  :  p.q  =  and  certain  aiixiliary  functions. 

U  pq 

Complete  details  are  given  in  [13],  where  proofs  of  the  properties  of 
^(•,h, t)  that  are  sted  in  the  following  lemma  caoi  also  be  found. 


Lemma  12.  There  exist  positive  constants  0  and  G^,  auid  continuous  func¬ 
tions  gQ  :  r"— and  :  {(h.t)  6  :  0  <  eh(l+||t||)  <  l}-^(0.oo) 

such  that 

Q 

i)  gpCt)  i  and 

i+iitir 

ii)  ||^(•,h,  t) -gQ(t)||^  g  i  Gj(h,t)h  for  each  h>0  and  t  e  r"^  that 
satisfy  0h(l+||t||)  <  1. 

A  consequence  of  (i)  is  that  we  ceuti  define  functions  u^  and  Uq(*:N), 
for  N  >  0,  in  H^(r'')  by 


Uq(x)  =  (2ir) 


-n/2 


f(t)gQ(t)e 


it*x 


dt 


R 


(23) 


Uq(x;N)  = 


(2ir) 


-n/2 


f(t)gQ(t)e 


it*x 


dt 


|tll<N 


whenever  f  €  L_(r”)  (cf.  [11]).  Note  that  u„(*;N)  =  u„  with  f  ®  f„.  We 
have  Uq  =  fg^  and  (Uq(*;N))  =  fj^gQ.  stnd  then  Parseval’s  equality  implies 

(24)  (|Uq-Uq(x;N)|1q  SG^IIf-f^ll^. 

The  integrals  that  define  u^  and  Uq(*; N)  are  to  be  interpreted  as  Lebesgue 
integrals  of  numerical -valued  functions.  However,  we  want  to  interpret 
Uq(*;N)  as  asi  integral  of  H^^iR*') -valued  functions.  In  the  appendix,  we 
show  that  the  integral  in  (23)  can  be  interpreted  in  both  ways,  in  aui 
unambiguous  and  consistent  manner.  Furthermore,  note  that  u^  is  the 
solution  of  the  constant  coefficient  differential  equation 
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(25) 


Above  we  have  used  expauislon  (22)  of  the  function  ♦(•,h,t).  Function 
♦(•,h,t)  Is  the  solution  of  the  elliptic  problem  (3)  with  h  auid  t  being 
parameters.  We  can  find  ♦(•,h,t)  numerically  e.g.  by  the  finite  element 
method  for  various  h  and  t  and  then  madce  various  expansions  or  approxima¬ 
tions  by  numerical  approaches.  The  a  posteriori  error  analysis  will  lead  also 
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to  an  assessment  of  the  accuracy  of  the  used  homogenization. 

Ue  can  use  various  approximations  of  4(*,h,t)  to  derive  the  homogeniza 
tion  approach.  Let  us  mention  a  few.  For  more,  see  [5]. 

a)  For  given  a  we  can  approximate  ♦(•,h, t)  so  that 

♦  (•.h.t)  a 

where  rational  functions  in  t.  Every  then  is  the  symbol 

of  the  homogenized  equation 

Lj(-)u  »  L2(*)u. 

b)  We  can  also  use  the  above  approximation  of  4(«,h,tl  as  an  "ainsatz" 

for  the  solution  and  derive  equations  for  via  finite  element  method 

auid  energy  minimization  principle.  This  approach  is  essentially  equivalent  t 
the  derivation  of  special  elements  for  "rough"  problems. 

c)  Using  other  approximations  we  can  obtain  pseudodifferential 
equations. 

d)  Any  homogenization  approach  can  and  has  to  be  understood  as  an 
approximation  of  ♦(•,h,t).  Hence  it  is  possible  to  aissess  the  accuracy  of  a 
homogenization  approach  by  comparing  it  with  '((•.h,  t). 

By  the  above  mentioned  approach  a  system  of  second  or  higher  order  el  lip 
tic  equations  can  be  derived.  This  system  will  usually  be  of  singularly  per¬ 
turbed  type  and  the  solution  will  show  boundary  larger  when  Q  is  bounded. 

We  can  use  various  "auisatzes"  of  various  types  in  different  areas  of  the 
domain,  especially  in  the  boundary  region  By  this  way  we  can  also  obtain 
reliable  solutions  in  the  neighborhood  of  the  boundary. 
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6. 


Define  w(x,t)  =  where  f  €  LgClR^)  and  is 

defined  in  Lemma  13.  According  to  (22),  Uq(*;N)  e  H^(IR^)  is  defined  by 


Uq(x;N)  = 


w(x, t )dt 


llt||<N 


as  an  integral  of  numerical-valued  functions.  Since 


ti — >  w(*,t)  is  strongly  measurable  in  and 

f  llw(*.t)||  dt  i  (2w)'^^G  |f(t)|dt  <  00, 

''||t||<N  ’  ^  ‘'l|t||<N 


it  follows  that 


W  = 


w( •, t)dt 


ItllSN 


can  be  defined  as  an  integral  of  H^^dR^'l-valued  functions.  We  want  to  prove 
Lemma  14.  Uq(*; N)  =  W  2ls  a  function  in  H^^(IR^). 

Proof.  It  suffices  to  show  that  xi — >UQ(x;N)e  and  xi — >W(x)e 

generate  the  same  generalized  function.  It  follows  from  (27)  that  a  sequence 
ti — >Wj^(*,t)  e  of*  simple  functions  can  be  chosen  so  that 


lim||w  ( ^t)  -  w(*,t)||  _  =0  a.e. 

k-»oo  ^ 


lim 
k-+oo  J 


|Wj^(*,t)  -w(«,t)||^  _^dt  =  0, 


|t||<N 


aind 


|Wk(*.t)||i  ^  -V 


Then  by  definition,  W  =  lim  f  w.  (•,t)dt  in  (R^). 

k-»co  •’lltllSN 
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W(x)e 


-v\x\ 


0(x)dx. 


□ 
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